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I briefly review the current status of chiral perturbation theory (CHPT) in the 
meson sector. Emphasis is given on the quest for higher precision. I discuss two 
examples: one where it is difficult to make a good prediction {Kj^ — ► 77^77), and 
where CHPT, even when pushed to higher orders, cannot yield an increase in 
accuracy. The second one is tttt scattering where a very sharp prediction can be 
made by combining CHPT at the two-loop level with dispersion relations. 



1 Introduction 

Chiral perturbation theory (CHPT) is the low-energy effective theory of the 
strong interactions. It yields a systematic expansion of any Green function in 
powers of quark masses and momenta, that automatically respects the Ward 
identities implied by th&j3hiral symmetry of the strong interactions. Its present 
form is due to Weinbergld and Gasser and Leutwylera. They showed the advan- 
tages of the effective-field-theory language over the direct implementation of 
the Ward identities as in the current-algebra framework. In particular, what 
was known before as a very difficult problem, the calculation of the unitarity 
corrections to a current-algebra result, was reduced to a routine loop calcula- 
tion in a well-defined framework. After the convenient tools of the effective 
field theory were made available, many processes have been calculated at the 
one-loop level: wherever possible, the comparison to the experimental data 
has shown a remarkable success of the method. 

In the early nineties several precision experiments on kaons, pions and 
etas, were approved, and brought the challenge of increasing the precision in 
CHPT calculations. Going beyond one loop became a necessity. The first 
complete two-loop calculation concerned the cross section for the two-photon 
annihilation into two neutral pions □. This beautiful and difficult calculation 
opened up the field of two-loop calculations in CHPT, that was very active 
during the nineties. In fact if we consider only the two-light-flavour sector, 
all the phenomenologically relevant calculations have already been completed, 
|SKhereas in the SU{3) framework, they have started only in more recent years 
□la. Moreover, in the purely strong sector, the Lagrangian at-prder and the 
complete divergence structure have been recently calculated 
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In an effective tfieory it is not obvious that a higher loop calculation may 
lead to higher precision: new constants come into play at each order, and since 
they are a priori unknown, they may induce an uncertainty which is as large 
as the effect that one is calculating. The quest for higher precision within an 
effective field theory is a particularly nontrivial one. In what follows I will 
discuss in some detail why it is nontrivial, and how, in some special cases, one 
can nevertheless succeed in making very accurate predictions. 

2 The Lagrangian of chiral perturbation theory 

At lowest order the Lagrangian of CHPT is remarkably simple^: 

£(2) ^ _(^^[/^M[/t + [/^t + ^[/t) + c{QUQU^) 

+ C2{\&D^U^D^U)+c:i{27L,lR) , (1) 

where 

D^JJ = d^,U - ir^U + iUl^ x = 2B(M + s + ip) 

M = diag(TO„,TOd,TOs) Q = e • diag(2/3, -1/3, -1/3) , (2) 

where U is the usual exponential of the matrix containing the octet of the 
Goldstone bosons, while r^, Z^, s and p are matrices of external fields. {A) 
stands for the trace of A. The simplicity of this Lagrangian is particularly 
remarkable in view of the variety of phenomena that it describes: the first two 
terms describe the strong interactions of the octet of the pseudoscalars between 
themselves and with external sources. The third, proportional to C describes 
the electromagnetic effects due to the exchanges of one virtual photon, while 
the remaining terms account for the effects generated by virtual exchanges of 
a, W OT a Z between quarks. 

The five constants which appear in the Lagrangian are not fixed by the 
symmetry but are easily determined in the comparison to the simplest of the 
observables: 

F^=F + 0{m) , M^o = B(m„ + ma) , 

n.r2 nr2 _ 0^2 C* .(1/2) 2 (3/2) 

Having fixed these constants, we can now calculate any other new observable, 
and make a prediction. For example the tttt scattering amplitude and the K13 

'^For simplicity's sake I do not write down explicitly the 27plet, nor the the CP-odd part of 
the weak Lagrangian. 
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form factors are given by: 



Ais, t, u) = ^al = = 0.16 , (4) 



{n-{p)\V^\K''{k))^{k+p)^f+{t) + {k-p),J^{t) , t^{k-pf 

f+{t) = l , f-{t)=Q . (5) 

2.1 The Lagrangian at next-to-leading order and beyond: strong sector 
In the strong sector the situation is the following: 

O(p^) ^ 10 new LEG, Gasser and Leutwyleri (84), Li , 
0(/) ^ 96 new LEG, Bijnens, G.G., Eckeri (99), d . 

For all the ten Li 's we have either a good phenomenological determination, or 
a sound theoretical estimate. In order to get a prediction at order p'^ in this 
sector one simply needs to do the one-loop calculation, remove the divergences 
with the help of the counterterms, and insert the numerical values of the finite 
parts from the most recent estimates. At this level the situation is almost as 
good as at order p"^. This is best seen with a couple of examples, like the 
TTTT scattering lengths and the K^s form factors, that w^b^ve seen before at 
leading order. The next-to-leading order corrections are 



0.16- [1.25] = 0.20 , (6) 



/+ (i) = /+ (0) (^1 + A+ ^ j ' ^+ ^ ^ ^'^^^ + A+ = 0.031 . (7) 

The formulae above neatly illustrate the beauty of the method: the S-wave 
scattering lengths and the Kf,j, form factors are expressible in terms of other 
observables (like the scalar and vector radii, and (r)y, respectively, the 

two Z3-wave scattering lengths, aj'^), plus other small corrections (indicated 
by the A symbols). The relation between these observables is exact up to 
order p"*, and will get corrections if we want to increase the precision of the 
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calculation. By inserting the most recent values for the observables on the 
right-hand sides, one gets the numerical predictions given above. We will 



discuss the situation for the tttt scattering lengths in detail below, ^s for the 



-ftTes form factors, the slope A+ is quite well known, and the PDGcl reports a 
number which is in excellent agreement with the one loop calculation above: 
A+ = 0.0286 ± 0.0022. 

If we move up to order the situation worsens considerably, and it be- 
comes impossible to write down expressions as nice as those in Eq. 
We have to give expressions in terms of the new low-energy constants C^'s, 
and since we know practically nothing about them, it is an obvious ques- 
tion whether these calculations represent an improvement at all in the final 
numerical predictions. Before discussing the matter in more detail in the fol- 
lowing sections, we can already anticipate an answer at the intuitive level: The 
low-energy constants parametrise our ignorance of the dynamics which is not 
explicitly present in the formalism. They contain contributions from all the 
physics above the scale of the pseudoscalar mesons - the most prominent effect 
is the one due to the lowest-lying resonances, as it has been explicitly shownB 
at order p*. Whenever a resonance or any other "high-energy" phenomenon is 
expected to influence substantially a low-energy quantity, one should expect 
a very important role of the low-energy constants, and a serious difficulty in 
making a precise prediction. On the other hand, observables dominated by the 
dynamics of the pseudoscalar, can be very well calculated within CHPT, and 
it makes good sense to push the calculations to order p^, as one can expect the 
contribution of the low-energy constants to be suppressed. Explicit examples 
of both situations are given below. 

2.2 The Lagrangian at next-to -leading order and beyond: weak and electro- 
magnetic sectors 

Both in the electromagnetic and weak sector, the situation becomes unpleasant 
already at next -to-leading order, where the number of constants immediately 
gets too large in comparison to the available experimental input: 

E.m. sector: 




0(p2) ^ C{QUQU^) , 



0(/) ^ 21 new LEG, Urech, Neufeld and Rupertsberger0 (93), K, 



Weak sector: 



0(p2) ^ c2{\eD^U^D^U) 

0(/) => 37 new LEG, Kambor, Missimer, Wyler, EckerQ (91), 
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In both cases, given the practical unfeasibility of a determination of all the 
available constants, the typical strategy has been to focus on a few selected 
processes where the same constants occur. One can then exploit the exper- 
imental information on some processes to make predictions in other. There 
are plenty of examples in the literature. A fairly comprehensive review of the 
situation in the weak sector (especially for the radiative decays), including a 
list of the new experimental measurements that are possible at a $ factory can 
be found in Ref. 113. 



3 Analysis oi K ^ tt^^i up to order p^. 

Assuming CP conservation the A{K]^ 7r'^77) is determined by two invariant 
amplitudes, A{s, v) and B(s, v), s = (qi + 92)^, = Pk • (^i ~ 92)7 where qi^2 
are the momenta of theptwo photons, and pk that of the kaon. At order p'^: 
B ^0. At order p-^JB: A^ 4/s(s - M^)F{s/M^) + . . ., and S = 0, where 
F{x) is a loop function generated by the tttt intermediate state in the s channel, 
that represents the dominant effect at this order, and the ellipsis stands for 
other less important contributions .|— Although the shape of the spectrum was 
nicely confirmed by the experiment til'llj, the branching ratio was a factor three 
too small: 

( (1.7 ±0.3) X 10-6 (NA31) 
BR = <^ (1.86 ± 0.60 ± 0.60) x lO^^ (E731) (8) 
[ 0.67 X 10-6 0(p4) , 

therefore requiring large 0{p^) corrections. The calculations at order p^ lil 
have considered only the (possibly dominant) corrections to the pion loops, 
and added to this a polynomial contribution; 

..l,„_M,.(^)...2Mi^,^^. 



s 



B^G 



where F(x) and G(x) also come from the tttt intermediate state in the s channel. 
To get into agreement with the experiment one needed to have a large and 
negative ay: BR = 0.83 x 10"^ ^j^h ay = and BR = 1.60 x 10"^ ^j^h ay = 
—0.9. Also for the spectrum, unitarity corrections alone were not sufficient 
(and actually worsened the comparison), while an improved agreement with 
the data is obtained only with ay ~ —0.9. 

The outcome of this 0{p^) analysis is therefore a clear need for a very large 
contribution from the polynomial part. Do we have a dynamical explanation 
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for such a large constant? I find it instructive to go back by seven vears and see 
how the situation appeared then: Cohen, Ecker and Pich in Ref. [3 described 
the situation as follows: "Several model estimates of ay have been made in the 
literature. A fair summary of those attempt is that we know neither the sign 
nor the magnitude of ay" More recently, however, D'Ambrosio and Portoles 
EZl have built a Vector Resonance Model that does indeed get the right sign 
and size for this constant: Uy^ ~ —0.72. This number is now in^mazing 
agreement with the one extracted from a fit to the most recent dataEfl. 

Although D'Ambrosio and Portoles estimate of ay was only a postdiction, 
it is reassuring to have an understanding of the size of this constant. In fact it 
is not the only case where one can find this relative size between the various 
contributions in the chiral expansion. A well-known analogous example in the 
strong sector is the vector form factor. Its Taylor expansion around s = 
is usually defined as Fv{s) = 1 + l/6(r2)yS + CyS^ + 0{s^). Cy vanishes at 
order p^, and can be predicted with no parameters at order p^. Nowadays it 
is known up to order p^llj: 
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M2 

'In^" + % 



' V2 
pi 



9607r2M2i^2 (i67r2F2)2 
(0.62 + 1.96 + 1.3 • 10"V^>2(Mp)) GeV"'' = 5.4 GeV 



(9) 



where the latter value is determined experimentally. Here also: i) the order 
parameter-free prediction fails badly; ii) there are large Oijfi) unitarity 
corrections; iii) but even larger Oijfi) polynomial contributions, coming from 
the well-known p resonance. 

4 TTTT scattering at order 

This is the "golden reaction" for Chiral Perturbation Theory: at threshold 
the naive exaansion parameter is M"^/! GeV^ ^ 0.02, and already a tree level 
calculation c2l should be rather accurate, in principle. However, this rule of 
thunib. is quite misleading here, as it is siupwn by the fact that both the one- 
loop (see Eq. (^)) and the two-loop E3 calculations produced substantial 
corrections. The violation of the rule of thumb has a well known origin, and 
is due to the presence of chiral logarithms L — MHIAttF-^Y In A/2//i2, which, 
for /i ~ 1 GeV change the expansion parameter by a factor four. If we look 
at the / = S'-wave scattering lengths, e.g., a large coefficient in front of the 
single (at one loop) and Hmihlp (^at two loops) chiral logarithms is the main 
source of the large correctionOa: 

™- + + 1^0.217 (10) 



327rF2 2 42 
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The final numbeiL.is the estimate given by the authors of the complete two- 
loop calculation c3, and was given without an estimate of the uncertainties, 
that required considerably more work. A better numerical prediction with an 
accurate estimate of the uncertainties is now available^, but before describing 
how this was obtained, it is useful to analyse in more detail the number in Eq. 

In particular it is important to check the role of the new (and largely 
unknown) Oijfi) constants. In this case their influence is rather small: to 
tain the number in Eq. ( [To| ) they were estimated using resonance saturation 
However, if one puts them to zero, the change in Oq is less than 1%. 
For the analytic part of the chiral expansion the rule of thumb works quite 
well. On the other hand it is more worrisome the effect of the uncertainties 
in the O(p^) constants, supposedly known quite well: changing the constants 
within their error bars changes the value of «[] typically by 0.01: The two-loop 
contribution can be completely overshadowed by the uncertainty coming from 
the 0(p*) constants - to improve the prediction one has to reduce drastically 
this uncertainty. This improvement in precision has been obtained by the use 
of a different method to determine the low-energy constants which we are now 
going to describe. 



4.1 Dispersive representation of the scattering amplitude 

As shown by RoyS, the fixed-t dispersion relations can be written in such a 
form that they express the tttt scattering amplitude in terms of the imaginary 
parts in the physical region of the s-channel. The resulting representation for 
j4(s, t, u) contains two subtraction constants, which may be identified with the 
scattering lengths Cq and Cq. Unitarity converts this representation intip a set 
of coupled integral equations, which was recently examined in detailed. The 
upshot of that analysis is that Oq and Oq are the essential low energy param- 
eters: Once these are known, the available experimental data determine the 
behaviour of the tttt scattering amplitude at low energies to within remarkably 
small uncertainties. 

The branch cut generated by the imaginary parts of the partial waves with 
t>2 starts manifesting itself only at 0{p^). Accordingly, we may expand the 
corresponding contributions to the dispersion integrals into a Taylor series of 
the momenta. The singularities due to the imaginary parts of the S- and P- 
waves, on the other hand, show up already at O(p^) - these cannot be replaced 
by a polynomial. We subtract the corresponding dispersion integrals as many 
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times as is needed in the chiral representation, for later convenience: 

W'(s)^ / ds' — r , (11 

where {eq, ci, £2} = {0, 1, 0}. Since all other contributions can be replaced by 
a polynomial, the phenomenological amplitude takes the form 

A(s, t, u)=16™g + {2al - bal) s + P(s, t, u) + 327r |i W^°(s) - ^W^is) 

+ - u)W\t) + |(s - + \W^{t) + + o(p8) _ 

We have explicitly displayed the contributions from the subtraction constants 
Oq and Cq. The term P[s,t,u) is a crossing symmetric polynomial 

P{s,t,u)^i+P2S+pj,s^ +Pi{t-uf +P5S^ +PQs[t~uf . (12) 

Its coefficients can be expressed in terms of integrals over the imaginary parts 
of the partial waves. In the following, the essential4ioint is that the coefficients 
Pi, . . . , P6 can be determined phenomenologicallyEj. 
In the chiral representation we have, similarly 

A{s, t, u)=C(s, t, u) + 327r { i C/°(s) + § (s - u) U\t) + | (s - t) U\u) 
+ lU^t) + ^U^u)~lU\s)}+0{p') , 

where C{s,t,u) is a crossing symmetric polynomial, 

C(S, t, u) = Cl + S C2 + C3 + (i — u)^ C4 + C5 + S (t — u)^ C6 , 

whose coefficients are given in terms of the low-energy constants @. Jlh£ 
functions U^{s), U^{s) and U^{s) describe the "unitarity corrections" E3t3 
associated with s-channel isospin / = 0,1,2, respectively, and are the chiral 
expansion (up to and including order p^) of the W^{s), Eq. ([III). 



4-2 Matching the two representations 

In their common domain of validity, the dispersive and the chiral represen- 
tations of the scattering amplitude agree, provided the parameters occurring 
therein are properly matched. Since the differences between the functions 
W-'{s) and C/^(s) are beyond the accuracy of the two-loop representation c3, 
the two descriptions agree if and only if the polynomial parts do, 

C{s, t, u) = 16^2 + (2^0 - 5al) s + P{s, t, u) + 0(/) . (13) 
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Since the main uncertainties in the coefficients of the polynomial P{s,t,u) 
arise from their sensitivity to the scattering lengths Cq, Oq, the above relations 
essentially determine the coefficients ci , . . . , Cg in terms of these two observ- 
ables. On the other hand, these coefficients in the chiral representation are not 
completely free, but are given in a power series in the pion mass. For example 
the combinations 

Cl = {C2 + 4M2(c3 - C4)} , C2 = ^{-Ci+ U4ic3 - C4)} , 

satisfy the following low-energy theorems 

Ci = l + ^(r^)5 + || + 0(a , C,^l + ^(r^)5 + ||£3-il|+0(a 

where ^ = (Af^/47rF^)^ = 0.01445. If we use the information on the scalar 
radius and a rough estimate for the constant £3 = 2.9 ± 2.4, then two of the 
Cj's become known, and we may use the set of equations to determine the 
two scattering lengths, and four of the c^'s. 

This leads to the following numbers for the chiral expansion of the / = 
S'-wave scattering length: 

a" = 0.197 0{p^) 0.2195 0{p^) 0.220 0(/) , 

which show a remarkably fast and good convergence. The difference of this 
procedure with respect to a direct evaluation of the scattering length is that 
here the chiral expansion is only applied in the unphysical region, where it 
converges best - the continuation to threshold to evaluate the scattering length 
is then provided by the dispersive representation, that, as shown inE^, at low 
energy does not show significant uncertainties. On this basis, we can safely 
conclude that yet higher orders will not modify the above predictions beyond 
the uncertainties present in this matching procedure. Once this are taken into 
account, the final result for the S'-wave scattering lengths is 

ag = 0.220 ± 0.005, 2ag - 5a2 = 0.663 ± 0.006 (14) 

The result (Q) relies on the standard picture, according to which the quark 
condensate represents the leading order parameter of the spontaneously broken 
symmetry. The scenario investigated in ref. l3 concerns the possibility that 
the Gell-Mann-Oakes-Renner formula fails, the second term in the expansion 

= 2Bm{ \ — ^^^3 + 0(^^)} being of the same numerical order of magni- 
tude or even larger than the first. Note that for this to happen, the value of IZ3I 
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Figure 1: The shaded region represervtei the intersection of the domains allowed by the old 
data and by the Olsson sunMule, seeE£| for details. The ellipse iniltflates the impact of the 
new, preliminary Ke^ dataEfl. The cross shows the result of Ref. Ej. The three diamonds 
illustrate the convergence of the chiral perturbation series at threshold. The one at the 
left corresponds to Weinberg's leading order formulae. The two black solid lines are the 
boundaries of the universal band £3 



must exceed the estimate we used by more than an order of magnitude. The 
constraints imposed on Oq, Aq by the available experimental information are 
shown in the figure. The ellipse represents the 68 % ccpfidence level contour 
obtained by ccmbining the new, preliminary Ke^ dataES with earlier experi- 
mental results c3. Concerning the value of Aq, the ellipse corresponds to the 
range 0.2 < aj] < 0.25. 

The figure shows that the values of Aq and aj] are stmngly correlated. The 
correlation also manifests itself in the Olsson sum ruleEiJ, which according to 
ref. Ea leads to 2ag - 5ag = 0.663±0.021 + 1.13Aa[]- l.OlAag, in perfect agree- 
ment with the result in eq. (0). This combination, however, is not sensitive 
to £3 - accurate experimental information in the threshold region is needed to 
perform a thorough test of the theoretical framewa£|k that undedies the calcu- 
lation. The forthcoming results from BrookhavencS, CERNc3'E2l and Frascati 
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will provide such a test. 



5 Conclusions 

We have reviewed recent developments in CHPT, emphasising in particular 
the problems that arise when pushing the chiral expansion to order . We 
have discussed a few examples that show different physical situations: in one 
case the low-energy constants at order have a dominant role, whereas in the 
other only a marginal one. In the former case CHPT cannot provide a solid 
prediction with a good control of the uncertainties, whereas in the latter case, it 
is possible to yield a very precise prediction. The case of tttt scattering has been 
described in some detail, and we have shown that by combining a dispersive 
representation and the chiral expansion one can bring the uncertainties down 
to the few percent level Eil. 
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